An Efficient Modified “ Walk On Spheres ” Algorithm for the Linearized Poisson - Boltzmann Equation
نویسنده
چکیده
A discrete random walk method on grids was proposed and used to solve the linearized Poisson-Boltzmann equation (LPBE) [1]. Here, we present a new and efficient grid-free random walk method. Based on a modified “Walk On Spheres” (WOS) algorithm [11] for the LPBE, this Monte Carlo algorithm uses a survival probability distribution function for the random walker in a continuous and free diffusion region. The new simulation method is illustrated by computing four analytically solvable problems. In all cases, excellent agreement is observed.
منابع مشابه
Efficient modified ‘‘walk on spheres’’ algorithm for the linearized Poisson–Bolzmann equation
A discrete random walk method on grids was proposed and used to solve the linearized Poisson– Boltzmann equation ~LPBE! @R. Ettelaie, J. Chem. Phys. 103, 3657 ~1995!#. Here, we present an efficient grid-free random walk method. Based on a modified ‘‘walk on spheres’’ algorithm @B. S. Elepov and G. A. Mihailov, Sov. Math. Dokl. 14, 1276 ~1973!# for the LPBE, this Monte Carlo algorithm uses a sur...
متن کاملMonte Carlo Methods for the Linearized Poisson-Boltzmann Equation
We review efficient grid-free random walk methods for solving boundary value problems for the linearized Poisson-Boltzmann equation (LPBE). First we introduce the “Walk On Spheres” (WOS) algorithm [1] for the LPBE. Based on this WOS algorithm, another, related, Monte Carlo algorithm is presented. This modified Monte Carlo method reinterprets the weights used in the original WOS algorithm as sur...
متن کاملRandom walk algorithms for estimating electrostatic properties of large molecules
The problem of computing the electrostatic energy of a large molecule is considered. It is reduced to solving the Poisson equation inside the domain and the linear Poisson-Boltzmann equation in the exterior, coupled by boundary conditions. Monte Carlo estimates for the potential point values, their derivatives, and the energy are constructed. These estimates are based on the walk on spheres and...
متن کاملA Spectral Monte Carlo Method for the Poisson Equation. a Spectral Monte Carlo Method for the Poisson Equation *
Using a sequential Monte Carlo algorithm, we compute a spectral approximation of the solution of the Poisson equation in dimension 1 and 2. The Feyman-Kac computation of the pointwise solution is achieved using either an integral representation or a modified walk on spheres method. The variances decrease geometrically with the number of steps. A global solution is obtained, accurate up to the i...
متن کاملEfficient and Accurate Higher-order Fast Multipole Boundary Element Method for Poisson Boltzmann Electrostatics
The Poisson-Boltzmann equation is a partial differential equation that describes the electrostatic behavior of molecules in ionic solutions. Significant efforts have been devoted to accurate and efficient computation for solving this equation. In this paper, we developed a boundary element framework based on the linear time fast multipole method for solving the linearized PoissonBoltzmann equat...
متن کامل